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ABSTRACT 

The instability of a quasi-Kepler flow in dissipative Taylor-Couette systems under the presence of an homogeneous axial magnetic field 
is considered with focus to the excitation of nonaxisymmetric modes and the resulting angular momentum transport. The excitation 
of nonaxisymmetric modes requires higher rotation rates than the excitation of the axisymmetric mode and this the more the higher 
the azimuthal mode number m. We find that the weak-field branch in the instability map of the nonaxisymmetric modes has always a 
positive slope (in opposition to the axisymmetric modes) so that for given magnetic field the modes with m> always have an upper 
limit of the supercritical Reynolds number. In order to excite a nonaxisymmetric mode at 1 AU in a Kepler disk a minimum field 
strength of about 1 Gauss is necessary. For weaker magnetic field the nonaxisymmetric modes decay. 

The angular momentum transport of the nonaxisymmetric modes is always positive and depends linearly on the Lundquist number of 
the background field. The molecular viscosity and the basic rotation rate do not influence the related of-parameter We did not find any 
indication that the MRI decays for small magnetic Prandtl number as found by use of shearing-box codes. At 1 AU in a Kepler disk 
and a field strength of about 1 Gauss the a proves to be (only) of order 0.005. 
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1. Introduction 



estimations we shall use the values 



The longstanding problem of the generation of turbulence in var- 
ious hydrodynamically stable situations has found a solution in 
recent years with the MHD shear flow instability, also called 
magnetorotational instability (MRI), in which the presence of 
a uniform axial magnetic field has a destabilizing eflTect on a 
diflTerentially rotating flow with the angular velocity decreasing 
outwards. 

According to the Rayleigh criterion an ideal flow is stable 
against axisymmetric perturbations whenever the specific angu- 
lar momentum increases outwards 



dR 



(R^^y > 0, 



(1) 



where (R, 0, z) are cylindrical coordinates, and Q is the angular 
velocity. Kepler rotation with D oc R~^ ^ is thus hydrodynami- 
cally stable against axisymmetric perturbations. It is not stable, 
however, under the presence of a uniform axial magnetic field 
(MRI). It is also not stable against nonaxisymmetric perturba- 
tions under the presence of a toroidal field (AMRI). In these 
cases the stability condition ([T]) simplifies to 



dry 

d^ 



>0. 



(2) 



For real fluids the instability condition looks more complicated. 
In order to be unstable the rotation must be fast enough and the 
magnetic fields must not be too weak or too strong. The lower 
magnetic limit is fixed by the electric conductivity of the plasma 
while the geometry of the disk determines the upper magnetic 
limit. As protoplanetary Kepler disks are so cool that even the 
action of MRI is in question ("dead zones"). In the present paper 
the lower limit is attacked with the calculations. For numerical 



r]^4' 10^^ cmVs, 



10"^^ g/cm^ 



(3) 



(see Brandenburg & Subramanian 2005) for which we have de- 
rived a minimum field strength of 0. 1 Gauss in order to excite un- 
stable axisymmetric modes in the protoplanetary disk (Riidiger 
& Kitchatinov 2005). This value is based on the calculation of 
the minimum magnetic Reynolds number 



Rm' = 



(4) 



for fixed values of the half-thickness H of the disk and the mag- 
netic Prandtl number 



Pm : 



(5) 



The vertical magnetic field is measured in terms of the 
Lundquist number 



(6) 



for which the Rm' is the absolute minimum for given Pm. Note 
that in the definitions of Rm' and the molecular viscosity does 
not appear. The viscosity only appears in the definition of the 
magnetic Prandtl number ([5]). 

The instability map Rm' = Rm^S') for the axisymmetric 
mode and for fixed Pm is very characteristic. For fixed Rm' ex- 
ceeding a minimum value of order 10 the instability exists be- 
tween an upper limit and a lower limit of which itself must 
exceed a minimum value of order unity (see Fig. 2 in RK05). 
In the following the curve formed by the upper limits is called 
the strong- field branch and the curve formed by the lower limits 
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is called the weak-field branch. Both the branches have oppo- 
site slopes. For given > 1 there is just one value of Rm' for 
marginal instability so that for all magnetic fields (above a lower 
limit) there is one rotation rate above which the axisymmetric 
MRI exists. For a certain . the Rm' takes its overall-minimum 

mm 

Rm^-^. For all Pm between 10"^ and 10^ the modes with the low- 
est Rm' (which are easiest to excite) are axisymmetric. 

The dependence of the curves of marginal instability on the 
magnetic Prandtl number Pm is only weak. For Pm < 1 there is 
no visible influence of Pm. Axisymmetric global MRI even ex- 
ists for very small Pm with one and the same magnetic Reynolds 
number (but, of course, the ordinary Reynolds number takes very 
high values). 

For Pm > 1 the behavior is diflTerent: both the critical Rm^-^^ 
and . grow with VPm. Hence, the scaling for Pm > 1 

mm o ' ' o 

switches to the parameters 



Ha: 



and 



Rm* = 



(7) 



(8) 



(Kitchatinov & Riidiger 2004). Both the critical rotation rate and 
magnetic field now run with ^Jvri instead with t] as it is true for 
Pm < 1. Both expressions are identical for Pm = 1. Again, the 
ratio of both quantities is free of the two dissipation parameters. 
The ratio of the linear rotation velocity and the Alfven velocity 
of the vertical field is called the magnetic Mach number Mm = 
U^/Va with 



Va = 



(9) 



At a distance of 1 AU and for a magnetic field of 1 Gauss the 
magnetic Mach number is of order 100. As the possible magnetic 
field should be smaller than 1 Gauss this value of Mm is certainly 
a minimum. 

The behavior of the nonaxisymmetric modes is not so well- 
known. We have learnt from the theory of the azimuthal mag- 
netorotational instability (AMRI) that too fast rotation always 
destroys the instability. The AMRI follows from the interplay of 
diflTerential rotation and an toroidal current-free magnetic field. 
It is basically nonaxisymmetric. For a given Lundquist number 
S there are two critical Reynolds numbers Rm. The instability 
is supercritical only between the two values. With other words, 
AMRI is excited by fast enough rotation but it is suppressed by 
too fast rotation. Both the weak-field branch and the strong-field 
branch have the same positive slope in the plane Rm over S. The 
reason for this expressive phenomenon is the smoothing influ- 
ence of diflTerential rotation on nonaxisymmetric magnetic per- 
turbations. 

The question arises whether also the nonaxisymmetric 
modes of standard MRI are finally supported by the diflTerential 
rotation. The answer has consequences for i) possible dynamo 
models but also ii) for the magnetic field amplitudes necessary 
for the excitation of nonaxisymmetric modes. Let us assume for 
a moment that the weak-field branch of the instability map ful- 
fills the condition 



Rm ^ Mmweak S 
then 



H 



> 



(10) 



(11) 



for the minimum seed field necessary for the excitation of 
the m = 1 mode. Mm^eak is the (positive) slope of the curve 
according to the calculations. The result does not depend on the 
actual value of the magnetic diflTusivity. 

The same is true for the strong-field branches shown in Fig. 
[2] One finds the corresponding magnetic Mach number Mmstrong 
much smaller than Mm^eak- For given basic rotation the flow is 
unstable against nonaxisymmetric perturbations if 

Mmstrong < Mm < Mmweak. (12) 

For Kepler disks the relation ( pT) can be read as 

< Mmweak', (13) 

with j3 = jioP/B^ as the plasma /3. Hence, the /3 must be rather 
small to excite nonaxisymmetric modes. The value of 400 used 
by Fromang et al. (2007) is so high that the corresponding mag- 
netic fields could be much too weak to excite nonaxisymmetric 
instability modes. Kitchatinov & Riidiger (2010) for a simplified 
Kepler disk find the rather moderate values Mmstrong - 2 and 
Mm^eak - 10. For such small values the unstable domain of the 
nonaxisymmetric modes is rather restricted. In the present paper 
we shall find that the numerical values of the critical magnetic 
Mach numbers for the weak-field branch for modes in infinite 
cylinders are much higher. 



2. The model 

The MRI has been found by considering the stability prob- 
lem of Taylor-Couette (TC) flows of magnetized ideal fluids 
by Velikhov (1959). Because of the simplicity of this geom- 
etry much work has been done to study the MRI for flu- 
ids between cylinders rotating with diff'erent angular velocities. 
Theory (Riidiger & Zhang 2001; Ji et al. 2001) and experiments 
(see Stefani et al. 2009) revealed the possibilities to realize the 
MRI even in the laboratory. 

In the following, two concentric cylinders (unbounded in z) 
of radii Rin and Rout are considered with the rotation rates flin and 
Qout- The rotation profile between the cylinders may mimic the 
Kepler rotation law, i.e. we shall fix Dout = 0.35Din. The value 
ensures that the cylinders with 7?out = rotate like planets. 
We call this radial rotation profile as quasi-Keplerian profile. 

The boundaries are assumed to be impenetrable, stress-free 
and perfect conducting, which are valid also for almost all liquid 
metal experiments in the laboratory. The cylinders are infinite in 
the axial direction. Fricke (1969) and Balbus & Hawley (1990) 
started to apply the instability to important astrophysical appli- 
cations. In order to be more close to the physics of accretion 
disks we also considered in the nonlinear simulations a pseudo- 
vacuum condition for the outer boundary. 

The governing equations are 



du 1 1 

— -h (mV)m = — WP -h vAu -h —curl BxB, 

ot p Jdo 



dB 



— = curl(M xB) + rjAB, 
ot 



and 

div u = div B 



0, 



(14) 



(15) 



(16) 



where u is the velocity, B the magnetic field, P the pressure, v 
the kinematic viscosity, and rj the magnetic diflfusivity. 
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The basic state is ur = = Br = B^p = 0, B^ = Bq = const. 



and 



RQ = aR-\- 



where a and b are constants defined by 

1 -y^a 



with 



T] = 



- Q /?2 



(17) 



(18) 



(19) 



The axial field amplitude is now measured by the Hartmann 
number 



Ha = 



BqD 



(20) 



D = Rout - ^in is used as the unit of length, v/D as the unit 
of velocity. The rotation Q is normalized with the inner rotation 
rate fli^. The magnetic Reynolds number Rm is defined as 



Rm 



(21) 



The Lundquist number S is defined by S = Ha • ^^m. Expressed 
with the Alfven frequency 11a it is 



S 
Rm' 



(22) 



All the calculations in this paper are done for a model with 
7/ = 0.5. This is the only geometry where the scales Ri^ and 
D are equal. As the magnetic Mach number the quantity Mm = 
Din/QA is used. 



3. Wave number and drift frequencies 

The equation s (p^ -p^ are linearized with respect to the back- 
ground state (|17|). The perturbed quantities are developed after 
azimuthal Fourier modes 



F = F{R)e 



i(kz+m(p+a)t) 



(23) 



The results are optimized in the wave number k. Only the solu- 
tions with those k are of interest for which the Reynolds num- 
bers take a minimum. All solutions with another k have higher 
Reynolds numbers. One can also show that the solutions with 
a certain positive k are always accompanied by a solution with 
-k but with the same Reynolds number and drift frequency (for 
given Ha and m). As the pitch angle of the resulting spirals is 
given by dz/dcp = -m/k it is clear that both the solutions have 
opposite pitch angles so that the solution is always a combina- 
tion of a left screw and a right screw. 

The proof is as follows. After elimination of both pressure 
fluctuations and the fluctuations of the vertical magnetic field, 
B'^, the linearized equations are 



duR ur im _ 

+ — + + ^kuz = 0, 



dR R 



R 



(24) 



1 dU(p 



dR^ R dR 
-i|mRe— — + CO 



U(j) 



R^ 



+ k' 



U(j)- 



2im _ 1 5 / 9 D , 



m 
T 



1 d Uz 1 I m 

R~d¥^¥~dR~\¥ 



+ k 



+ --Ha 

k 



R 
IOBr 



Br 



R dR R^ 



+ ^Ha2 

k 



m 



(25) 



1 



1 du^ 



dR^ R dR^ R^ dR 

In? A D \duz . d I ^ , 
+ — - 1 mRe- + ^ - - imRe- - 



-Ha 



O^Br IOBr Br 

- k Br+ 



dR^ R dR 
im dB(p im 



R dR R^ 

^2 



ik 



R^ 

d^UR 1 dUR Ur 

R~dR ~ W 



-[k^ + ^^^ur - ki^Rt^ + cd^ur - 



km Q 
-2—rU0- 2ikRQ—-U0 = 0, 

R iZin 



(26) 



2im 



B(p - iPm I mRe 



d^B. 1 OB. 



dR^ R dR 



B(i) 



m 



-\- clj\Br-\- iku, 



%R 



(27) 



2im , I a 
-h — —Br - iPm mRe—— cl) 1 5^ ikU(p -\- 

^FmReR—^BR = 



(28) 



(Shalybkov et al. 2002). The eigenfrequency oj is here normal- 
ized with the viscosity frequency v/7?q. One finds the system 
as invariant against the simultaneous transformation k -k, 
Uz -Uz, Br -Br and B^p -B^. Hence, if a solution is 
known for a certain k, then always a modified solution exists for 
-k. The standard MRI in cylindric geometry always produces 
the same number of left and right screw^ The total helicity (ki- 
netic and magnetic) is thus vanishing (Fig.[T]). 

Conjugating the linearized equation system ([24]-[28]) leads to 
the finding that the system is invariant against the simultaneous 
transformations m -m, k ^ -k and %{(jS) -%{(jS). Hence, 
both the pattern drift dcpldt = -%(oj)/m of both solutions and 
also the pitch 



dz 



m 
T 



(29) 



are equal so that both the solutions are identical. It is thus enough 
to consider the solutions for positive m. 



^ The same is true for the instability of toroidal fields between the 
cylinders without and with electric current, see Riidiger et al. (2011). 
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Fig. 1. The simultaneous existence of the two modes with k and 
-k in the nonhnear regime. The whole pattern drifts in the direc- 
tion of rotation. Rm = 160, S = 13, Pm = 1. 



The unstable Taylor-Couette flow forms axisymmetric or 
nonaxisymmetric vortices. With our normalization the vertical 
extent dz of a vortex is given by 



5z 



^out ^ii 




hence for fj = 0.5 



71 
k' 



(30) 



(31) 



For ^ ^ TT the cells have the same vertical extent as they have in 
radius and for ^ » tt the cells are very flat. Generally, the vor- 
tices of the axisymmetric modes become more and more elon- 
gated in the vertical direction, k ^n. 

The drift velocity "Rico) which always proved to be negative, 
i.e. the pattern drifts are in direction of the rotation (eastward). 
It is 



(32) 



so that the eastern drift period in units of the rotation period is 
m/'R(oj). Note that we are working in the resting laboratory sys- 
tem. 



4. Solutions for Pm=1 

We start with the standard case Pm = 1 where no diflference exist 
between the both Reynolds numbers and also no diflference exists 
between the Hartmann number and the Lundquist number. Most 
numerical simulations concern this case. The critical Reynolds 
numbers for excitation of the modes with m = 0, m = 1 and 
m = 2m quasi-Kepler flows are given in Fig. [2] 

For the given Lundquist numbers the Reynolds numbers are 
minimized by variation of the axial wave number k. Note the ex- 
istence of an absolute minimum Rcmin of the Reynolds number. 
It is smaller for the axisymmetric mode than for the nonaxisym- 
metric modes. The nonaxisymmetric modes need faster rotation 
for their excitation. There is a basic diflference, however, between 



1000 



E 
a: 








10 



20 



30 



40 



Fig. 2. The instability map of the modes with m = 0, m = 1 
and m = 2 in quasi-Kepler flows (jiq = 0.35). For given Rm the 
flow in unstable for Lundquist numbers S between the values at 
the weak-field branch and the strong-field branch. Star symbols 
mark the series of nonlinear simulations discussed below. Pm = 
1. 



the axisymmetric and the nonaxisymmetric modes. For m = 
and S ^ 1 exists always one critical Reynolds number above 
which the MRI is excited for all larger Rm. The absolute mini- 
mum value for S is of order unity. On the other hand, the critical 
Lundquist values for the nonaxisymmetric modes with m > 
behave completely diff'erent. For S > Smin where Smin - 1 is 
the smallest possible Lundquist number there are always two 
critical Reynolds numbers between them the nonaxisymmetric 
MRI modes can only exist (Fig. [2]). Hence, the nonaxisymmetric 
mode cannot survive if the rotation is too fast. The diff'erential 
rotation excites the MRI but - if too strong - it suppresses its 
nonaxisymmetric parts. 

The rotational quenching of the nonaxisymmetric parts of the 
MRI should have serious consequences for the dynamo problem. 
After the theorem by Cowling all dynamo models need nonax- 
isymmetric parts of the magnetic field. 

The m = 1 mode in Fig. |2] needs seed fields which are the 
higher the faster the rotation is. As an estimation one finds 



3 • 10"-^ a, 



(33) 



(Mm^eak - 320). If we take the linear velocity of the Kepler disk 
at 1 AU as more than 30 km/s (solar systern) and the density at 
this place as 10"^^ g/cm^, then the relation (33 ) means 



Bo > I Gauss, 



(34) 



which is a very high value at the distance. The value (34) is 
needed to start any MRI-dynamo. It exceeds the above men- 
tioned minimum value of 0.1 Gauss for the excitation of axisym- 
metric modes by one order of magnitude. The rather high value 
of this field strength suggests that indeed the weak-field branch 
of the MRI bifurcation map is the branch of astrophysical rele- 
vance rather than the strong-field branch. For the latter one finds 
Qa - 0.3 Din representing magnetic fields which are by a factor 
of hundred (!) stronger than the fields at the weak-field branch. It 
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is hard to imagine that such fields are available for the formation 
of protoplanetary disks. 



201 ' ' ■ ' 1 0.5 




10 20 30 40 50 10 20 30 40 50 

S S 



Fig. 3. The same as in Fig.[2|but for the normalized wave number 
(left) and the drift frequency (right). 

The results for the wave numbers (for which the Reynolds 
numbers for fixed S are minimum) are plotted in Fig. [3] (l eft). 
The dotted line gives k = n for which the cells after ( |31| ) are 
nearly spherical. For ^ < tt the cells tend to be prolate while they 
are oblate for k> n. This is, of course, only true if the radial cell 
size is of the order of the distance between the cylinders which 
must be checked separately. 

There are surprising diff'erences in Fig. [3] (left) for axisym- 
metric and nonaxisymmetric modes. All the axisymmetric rolls 
are prolate as their axial size exceeds the radial scale. This is true 
for all values of the magnetic field. The result is not surprising as 
magnetic fields always increase the correlation length along the 
field, hence for strong fields we have k oc 1/S. That for m = the 
wave numbers k for given Rm are small at both the weak-field 
(left) branch and the strong-field (right) branch of the marginal 
instability curve does not mean that they are small also between. 
Kitchatinov & Riidiger (2004) find with a local approximation 
close to the left branch the relation k ^ Va/tj so that the wave 
number runs with S. One must thus expect that between the left 
branch and the right branch there is a maximum of k so that 
there also the axisymmetric modes have short axial scales - as 
we shall demonstrate by Fig. [8] (top, right). Figure |4] gives for 
Rm = 500 the mentioned maximum of the axial wave number 
together with the growth rate. One finds that the axisymmetric 
channel mode forms thin rolls only close to the left branch but 
the rolls are thick for stronger fields, i.e values of Mm > 1. 

For m = 1 only the cells with large S are almost rolls. The 
wave numbers for both branches are rather diff'erent. They are 
large for the weak-field branch and they are small for the strong- 
field branch. The cells for the marginal modes along the weak- 
field branch are thus rather flat. This is only true, however, if the 
radial eigenfunctions do smoothly cover the distance between 
the cylinders. 

The drift of the nonaxisymmetric modes is eastward typi- 
cally 50% of the rotation frequency of the inner cylinder, i.e. the 
MRI pattern drift is westward with respect to the rotating sys- 
tem. As the outer cylinder only rotates with 35% of the rotation 
frequency, the field pattern rotates by 42% faster than the outer 
cylinder. The corotation radius of the MRI pattern is located be- 
tween both the cylinders close to the middle between the cylin- 
ders. The same is true, at least qualitatively, for the pattern with 
m = 2. A test calculation with ^ = 0.8 leads to a pattern drift 
-%(co)/m ^ 0.85 which also yields a corotation radius just in 
the middle of the two cylinders. 




10 20 30 40 50 
S 

Fig. 4. The dependence on the magnetic field of the wave number 
of the axisymmetric mode and the growth rate (dashed) for given 
Reynolds number Rm = 500. m = 0, Pm = 1 . 

5. Solutions for Pm=0.01 

The results for smaller magnetic Prandtl numbers (Pm = 0.01) 
are given in the plots [5] and |6] In the weak-field limit the diflTer- 
ences to Pm = 1 are rather small. The rotational quenching of 
the nonaxisymmetric modes seems to be slightly weaker than for 
larger Pm (Mm^eak - 400). 

Also the behavior of the wave number and the azimuthal drift 
is rather obvious. The cells in the weak-field limit are flat while 
they are highly aligned along the rotation axis in the strong-field 
limit (Fig. [6j left). The axisymmetric cells are longer than the 
nonaxisymmetric cells. 




10 20 30 40 

S 



Fig. 5. The same as in Fig. [2| but for Pm = 0.01. At the strong- 
field branch for S > 20 the nonaxisymmetric mode becomes the 
mode with the lowest magnetic Reynolds number. 

More striking is the phenomenon for strong fields that the 
marginal- stability curves for m = and m = 1 are crossing for 
S ^ 18. The mode with the global-minimum Reynolds num- 
ber is always axisymmetric but for stronger magnetic fields the 
critical Reynolds numbers for m = 1 are smaller than those for 
m = 0. For such fields and for increasing Reynolds numbers MRI 
sets in as a nonaxisymmetric flow pattern. The nonaxisymmetric 
structure is lost, however, for too fast rotation when the mag- 
netic Reynolds number reaches the upper value of the marginal 
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Fig. 6. The normalized wave number (left) and the drift fre- 
quency (right) for Pm = 0.01. 



instability of the m = I mode and the solution becomes axisym- 
metric again. We have found this sort of mode crossing both for 
different geometries (Kitchatinov & Riidiger 1997) and also with 
different codes (Shalybkov et al. 2002). So far this phenomenon 
is only known for MHD flows between conducting cylinders. 

6. The magnetic Prandtl number dependence 

There are simple scaling laws for the marginal instability. They 
appear if also the instability domains for large magnetic Prandtl 
numbers are computed. The results for Pm = 10 are plotted with 
the variable Rm* and Ha (Fig. [7] right). Both the Figs. [7] demon- 
strate the small differences of the curves for small Pm if scaled 
with Rm and S, while the same is true for the curves for large Pm 
if scaled with Rm* and Ha. In both cases the ratio of the char- 
acteristic numbers is the magnetic Mach number which is free 
of any diflTusivity values. We find a somewhat weaker rotational 
quenching 



Da ^ 3 • 10""^ Qi, 



(35) 



of the nonaxisymmetric modes for Pm > 1 compared with the 
relation ( [33] ) which holds for Pm < 1 . It is the general behavior 
of the weak-field edges of the instability domains for Pm < 1 to 
show a positive slope. However, after inspection of Fig. [7] (right) 
for Pm = 10 the slope is very large and it is unclear whether 
it is still positive. In any case the behavior of the weak-field 
limit for large Prandtl numbers slightly differs from that of the 
curves for smaller Pm. The differences provided by the strong- 
field branches are very small for the two regimes. 




Fig. 7. The critical Reynolds numbers for the nonaxisymmetric 
modes with m = 1 for small (left) and large (right) magnetic 
Prandtl numbers. 



7. Nonlinear simulations 

With the 3D spectral MHD code for cylindric geometry de- 
scribed by Gellert et al. (2007) also nonlinear simulations of 
global MRI are possible. The code works with M Fourier modes 
in the azimuthal direction, a typical value used in the calculations 
is M = 16 for the rather narrow spectrum of excited modes. 



7.1. The instability pattern 

The simulations concerning the instability pattern are related to 
the map of marginal instability for Pm = 1 (Fig. [2]) for a fixed 
Lundquist number S = 13. Three examples are given. The first 
one works with relatively slow rotation while the second con- 
tainer rotates faster or much faster. The marginal instability ap- 
pears for a minimum magnetic Mach number of about 6. The 
first model lies below the instability domain of m = 1, the sec- 
ond one inside the m = 1 domain and the third model lies out- 
side. The magnetic Mach numbers for the considered cases are 
Mm = 7, Mm = 46 and Mm = 96. 

The nonlinear calculations are important as our curves for 
marginal instability of nonaxisymmetric modes only concern the 
stability behavior of the Kepler flow under the presence of the 
axial field. It is also possible that the rolls of the m = solu- 
tion become unstable against disturbances with m> leading to 
secondary instabilities with nonaxisymmetric patterns. 

Figure [8] shows the results. For given magnetic field ampli- 
tude (S = 13) simulations with Rm = 88, Rm = 600 and one 
with Rm = 1250 are presented. Only the value Rm = 600 lies 
in the instability map (Fig. [2]) between the lower and the up- 
per limit for nonaxisymmetric instability and only in this case a 
nonaxisymmetric (drifting) magnetic pattern results ([8j middle). 
From the beginning on, the mode with m = grows fastest and 
also the mode with m = 1 grows continuously. The other modes 
come much later so that only after 30 orbits the complete pattern 
occurs. 

The opposite is true for faster rotation (Fig. [8] right). Here 
only the mode with m = grows while all the nonaxisymmetric 
modes decay (Fig.[8j bottom right). The resulting magnetic pat- 
tern is axisymmetric despite of the high value of the Reynolds 
number. It is thus shown that indeed the linear approximation 
provides the real instability behavior. For the considered rota- 
tion rates a nonaxisymmetric instability of the magnetic m = 
pattern does not appear in the numerical simulations. 

Note that - as predicted by the linear results for the ax- 
ial wave numbers (Fig. [3] left) - the axisymmetric cells for 
Reynolds numbers close to the instability limit are really prolate 
(Fig. [8] left). They do become oblate for the fast-rotation case 
(Fig.[8| right) despite of the action of the Taylor-Proudman the- 
orem. The nonaxisymmetric modes form a spectrum of modes 
(see Fig. [8] middle). Only models of this kind are used in the 
next Section to compute the outward angular momentum trans- 
port of the MRI. Only for such models the equations form a non- 
linear mixture of many azimuthal modes close to the transition 
of turbulence. Obviously, too fast rotation destroys the resulting 
mixture (Fig. [8] bottom right). 

The amplitudes of the MRI-induced magnetic fluctuations 
are also of importance. They easily exceed the strength of the 
axial background field but only in the domain where the nonax- 
isymmetric modes are excited (see Fig. [8]). On the other hand, 
the amplitudes of the toroidal field components grow with grow- 
ing Reynolds number. For Re = 600 the maximum B^p exceeds 
the axial background field by one order of magnitude. 
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Fig. 8. Top: the radial component of the magnetic field for S = 13 and for slow rotation (Rm = 88, left), medium rotation (Rm = 600, 
middle) and fast rotation (Rm = 1250, right). For slow and for fast rotation the nonlinear instability pattern is purely axisymmetric 
and for medium rotation it is nonaxisymmetric. Bottom: the same for the energies in the Fourier modes of the magnetic field. 
Pm= 1. 



7.2. The angular momentum transport 

The total angular momentum transport in radial direction is 

which in the usual manner can be written as 



ClK 



(37) 



with D as the gap width of the container. Hence, the MRI-ck as 
the normalized angular momentum transport can be computed 
with the definition 



(38) 



Note that this definition diff'ers from the usual one unless D ^ H, 
which values are indeed of the same order for thick disks. In 
the first step we computed the quantity ([38]) by averaging only 
over the azimuth (Fig. [9]). One learns that the angular momentum 
transport is positive everywhere with a weak indication of the 
cell structure. There are no areas in the computational domain 
with negative angular momentum transport. Generally, after our 
experiences the Maxwell part in the relation ( 36 ) dominates the 
Reynolds term. 

In order to overcome the axial inhomogeneity we continue 
to average over the z coordinate so that the resulting expression 
remains only a function of R. The resulting profile shows a char- 
acteristic maximum as it must vanish at the boundaries because 
of the boundary conditions. The question is how this maximum 
is related to the mean pressure in the computational domain. 
The pressure in the container vanishes at the inner cylinder and 




Fig. 9. The angular momentum transport parameter a if averaged 
only over the azimuth is everywhere positive and only slightly 
reflecting a cell structure. Rm = 400, S = 30, Pm = 1. 

monotonously grows towards the outer cylinder. The radially av- 
eraged pressure for Kepler rotation between the bounding cylin- 
ders is of order pQ^D^ so that indeed the relation (37) can be 
read as 



pTr^ = ass P- 



(39) 



This is the standard relation of the accretion theory which - as 
we now know - only holds in the model after averaging over the 
entire cylinder. 
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In the next step, therefore, the averaging procedure concerns 
the full container. The following results for the ass parameter do 
only concern to this model. They are given in Fig. [T0| and they 
can be represented by the linear relation 



Qfss = 4.8 • 10"^ S, 



(40) 



valid for all the considered Reynolds numbers and magnetic 
Prandtl numbers (Pm < 1) given in Fig. [2] with star symbols. 
There is no dependence of the ass on the rotation rate except the 
chosen value of the magnetic Reynolds number lies too close 
to the boundaries of the instability map (Fig. [2]). For two exam- 
ples for Pm = 1 (green diamonds in Fig.[TO]) the outer boundary 
condition has been changed from perfect conductor to pseudo- 
vacuum. The numbers do not show a remarkable influence of the 
boundary conditions on the resulting values of ass • 




our ass as a realization of the 13 viscosity in the sense of Duschl 
et al. (2000) and Hure et al. 2001). It is obvious that the magnetic 
field amplitude must not much smaller than about one Gauss in 
order to get ass values of order 0.01 or even larger. 

The angular momentum transport shown in Fig.[TO|is only 
due to the nonaxisymmetric modes with m > 0. Only these 
modes have been defined as the 'fluctuations' in the definitions 
of the random functions u and B and the well-defined averag- 
ing procedure is considered as the integration over cp. We have 
shown in Fig. [8] that only these modes in our simulations are 
close to develop turbulence. 

The question remains whether also the axisymmetric modes 
with m = contribute to the angular momentum transport. They 
can be defined as fluctuations only by their small-scale spatial 
vertical variations characterized by the vertical wave number k. 
Then the basic rotation can only be defined by (w^) after aver- 
aging over z. The fluctuations of the flow and the field defined 
in this way do indeed transport angular momentum. If also the 
axisymmetric modes are used for the calculation of the angular 
momentum then the resulting values always lie above the curves 
for the nonaxisymmetric modes (Figs. [TT] and [12]). For higher 
Reynolds numbers the extra value by the axisymmetric modes 
diminishes so that the lines in the Figs. [TT] and [12] are approach- 
ing. In this picture the relation (40 ) remains true if the Reynolds 
numbers are large enough compared to the Lundquist number 
S, i.e. in the regime of high magnetic Mach number. With such 
large values of the magnetic Mach number as discussed above 
the contributions of the axisymmetric modes to the viscosity-of 
indeed remain small. 



500 



1 000 1 500 2000 
REYNOLDS NUMBER 



2500 



x10"' 




30 40 
LUNDQUIST NUMBER 

Fig. 10. The normalized angular momentum transport ass (aver- 
aged over the entire cylinder) in its dependence on the Reynolds 
number (top) and the Lundquist number S (bottom). Influences 
of both the viscosity and the basic rotation on the angular mo- 
mentum transport parameter ofss do not appear. Dots: Pm = 1.0, 
circles: Pm = 0.5, stars: Pm = 0.2, square: Pm = 0.1. Green di- 
amonds: Pm =1.0 and pseudo- vacuum as outer boundary con- 
dition. 



From Eq. (40) one finds with the numbers (3) that 
B D 



ass ^ 0.005 



IGauss lAU 



(41) 



Hence, the numerical value of the MRI-(2 linearly depends on the 
amplitude of the magnetic field and/or the size of the disk or the 
torus. According to the definition ([37]) one can also understand 
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Fig. 11. The normalized angular momentum transport ass (aver- 
aged over the entire cylinder) in its dependence on the Reynolds 
number. Open circles represent the values if the contribution of 
the channel mode (m = 0) is added. S = 13, Pm = 1. 



It is shown in Fig. 1 1 that the angular momentum transport 
by the nonaxisymmetric modes stops for Rm ^ 1200. For larger 
Reynolds numbers the axisymmetric rolls shown in Fig. [8] (top, 
right) do alone produce the viscosity-(2 with almost the same 
value of about 10"^. The given Reynolds numbers reach the nu- 
merical limits of the code so that the angular momentum trans- 



8 



M. Gellert, G. Riidiger & M. Schultz: Standard MRI in quasi-Kepler Taylor-Couette flows 




CO 2 



500 


1000 


1500 








xlO"^ 










-•- averaged - 






-e-incl. m=0 




500 


1000 


1500 



REYNOLDS NUMBER 
Fig. 12. The same as in Fig.[ll]but for S = 30, Pm = 1. 

port by the m = mode for very large Rm remains an open 
question as well as its stability towards larger Rm. 

Our results do not confirm the finding of Lesur & Longaretti 
(2007) and Longaretti & Lesur (2010) who report strong depen- 
dencies of the normalized angular momentum transport on the 
microscopic viscosity and/or the basic rotation rate. In particu- 
lar we do not see a decay of the angular momentum transport 
of the MRI for small magnetic Prandtl number in contrast to 
the conclusions of Fromang et al. (2007). As known from the 
bifurcation maps of the standard MRI the actual value of the 
magnetic Prandtl number is not very important if only the basic 
rotation and the magnetic background field are normalized with- 
out use of the microscopic viscosity. Of course, one can split the 
Lundquist number in accordance to 



Rm 



(42) 



with the plasma p = (flD ^|^^|Bzf'. But in this case the depen- 
dence of a on Rm and p such as described by Longaretti & Lesur 
(2010) is merely an artificial dependence. After our findings the 
Qfss only depends on the linear product of the magnetic field, the 
electric conductivity and the disk size. The total value of ass 
does not strongly exceed the limit of 10"^. In comparison with 
the results of Longaretti & Lesur the global simulations lead to 
values smaller by one order of magnitude. 

8. Discussion 

The standard MRI is the instability of diff'erential rotation un- 
der the presence of an axial magnetic field. Axisymmetric and 
nonaxisymmetric modes are excited for supercritical but not too 
strong fields. One must expect, however, that the nonaxisym- 
metric modes are quenched by too strong diff'erential rotation. 
On the other hand, the nonaxisymmetric modes are important 
for any form of hydromagnetic dynamo and/or the evolution of 
small-scale turbulence and its angular momentum transport. It is 
this rotational quenching of the nonaxisymmetric modes excited 
close to the weak-field branch and its consequences for the an- 
gular momentum transport which is probed in the present paper. 



We have shown by considering a dissipative cylindric TC 
flow that the MRI only forms nonaxisymmetric modes in a spe- 
cial domain in the Reynolds number-Lundquist number map. 
While for large Rm the minimum Lundquist number Smin - 1 
for the excitation of axisymmetric modes does not depend on 
Rm the minimum magnetic field for the excitation of nonax- 
isymmetric modes grows for growing Rm. Both the weak-field 
limit and the strong-field limit of MRI have a positive slope for 
nonaxisymmetric modes (see Kitchatinov & Riidiger 2010). The 
inverse magnetic Mach number I^A/^in must exceed the value 
^ 3 • 10"^ (for Pm = 1) in order to get nonaxisymmetric modes 
excited. This quenching eff'ect is very similar for Pm < 1 but it 
is weaker for Pm > 1 (see Eq. (35 )). Fluids with Pm > 1 seem 
to be favored, therefore, in numerical simulations. 

For field amplitudes smaller than the given ones the standard 
MRI in cylinders is only formed by axisymmetric rolls which re- 
mained stable in our simulations with Rm up to 10^. The domain 
of turbulence, therefore, proved to be (much) smaller than the 
entire domain of MRI. The most important numerical result of 
the simulations is that for increasing magnetic Reynolds number 
the MRI develops from axisymmetric rolls to nonaxisymmetric 
'turbulence' and returns to axisymmetric rolls. We did not find 
sofar a second bifurcation of the roll instability to more chaotic 
patterns (Fig. [8]). 

The dependence of the wave numbers on the field strength 
is in accordance to the expectation that the axial scales grow for 
growing field. Figure [3] shows for m > how large values of S 
correspond to small values of k as it follows from the Taylor- 
Proudman theorem. The same is true for m = but with a spe- 
cialty: close to the weak-field branch the axial wave length sinks 
down to a minimum for increasing field amplitudes due to the 
influence of the dissipation. Beginning from this minimum for 
increasing field amplitudes the wave lengths grow in accordance 
to the Taylor-Proudman theorem (see Fig.|4]). 

The angular momentum transport by the MRI has been com- 
puted with a nonlinear MHD spectral code which, of course, is 
able to reproduce all the above given results of the linear the- 
ory. In a first step only the angular momentum transport by the 
nonaxisymmetric modes (which in the model are the only one 
to be responsible for turbulence) has been computed. Expressed 
by the MRI-a (38) the results are of a striking simplicity. We 
did not find an influence of the numerical values of the mag- 
netic Prandtl number and/or the magnetic Reynolds number on 
the resulting a^ss (Fig. 10 ). The Lundquist number gives the only 
influence represented by the linear relation ( [4Q| . The numerical 
value of the ass taken at (say) S ^ 100 (1 Gauss for proto- 
planetary disks at R= 1 AU!) is only of order 10"^. This value 
strongly reminds on the numerical results of Brandenburg et al. 
(1996) obtained with shearing box simulations. Higher values of 
ass require stronger fields. Nonlinear relations between a^ss and 
the given axial field component S have not been confirmed here. 

The inclusion of the angular momentum transport by the axi- 
symmetric 'channel' modes gives only a slight modification of 
this picture. The maximum amplification of the angular momen- 
tum transport by the channel mode is characterized by a factor 
of two, which for slow rotation gives a stronger increase of the 
Qfss than for fast rotation (Figs. [TT]and[T2]). The total angular 
momentum transport expressed by Ofss due to axisymmetric and 
nonaxisymmetric modes finally obtains a weak decrease for in- 
creasing Reynolds number. There is even indication that for large 
Reynolds number when the nonaxisymmetric modes eventually 
become stable the transport by the axisymmetric modes yields 
nearly the ass as the transport by the nonaxisymmetric modes 
for medium Reynolds numbers does. 
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